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A new test apparatus is constructed and is applied to investigate static and 
dynamic characteristics of annular seals for turbopumps. The fluid forces acting on 
the seals are measured for various parameters such as the preswirl velocity, the 

pressure difference between the inlet and outlet of the seal, the whirling 
amplitude, and the ratio of whirling speed to spinning speed of the rotor. 
Influence of these parameters on the static and dynamic characteristics is 
investigated from the experimental results. As a result, preswirl affects the 

dynamic characteristics strongly. Especially, the preswirl opposing the rotating 

direction has a stabilizing role on the rotor system. 

1. Introduction 

A turbomachine such as a pump tends to be operated at high speed and high 
pressure because of severe performance requirements. But with increasing speed and 
pressure, the fluid force causes an instability of the rotating machine because of 
the noncontacting seal, such as the balance piston and wear ring. So, it is 
necessary to make the dynamic behavior clear and to supply accurate data in order to 
predict and prevent the unstable vibration, and still more necessary to design a 
stable rotating machine. 

Many theoretical and experimental analyses on the characteristics of 

noncontacting seals are reported. As an experimental research on annular seals, 
Childs et al. (Refs. 2,3) investigated fluid forces acting on the seals and their 
characteristics with a test apparatus which has an eccentric rotor. Nordmann (Ref. 
4) used an impulse force type test apparatus to carry out an experiment on the 

characteristics of annular seals. Kaneko et al. (Ref. 7), Kanki et al. (Ref. 8), 
and Hori et al. (Ref. 9) also reported their test results on the seals. All these 
results are important for the research on annular seals. However in these 
experimental apparatus, the casing inertia force of the test rig is larger than the 
flow-induced force due to the seal, so the signal/noise ratio is not good. Childs 
et al. investigated the characteristics under the condition of rotor whirling 
motion, but limited by the test apparatus, he only investigated the characteristics 
in the synchronous motion of the rotor. Recently, Kanemori and Iwatsubo (Ref. 10) 
used a test apparatus in which the rotor gives spinning and whirling motion 
independently to test a long annular seal of a submerged pump for various whirling 
speeds. But, they did not test for the preswirl velocity, which strongly affects 
the stability of a rotor system, because their interest is limited to the dynamic 
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characteristics of the submerged pump. The dynamic characteristics of annular seals 
have not been investigated yet in various operating conditions. 

This paper shows the experimental results obtained by an experimental apparatus 
which is newly designed to obtain a high S/N ratio. In the test the fluid force 
acting on the seals is measured and the stiffness, damping, and inertia coefficients 
are obtained. Influence of the parameters such as the whirling amplitude pressure 
difference between the inlet and outlet of the seal, preswirl, rotating and whirling 
speeds, and their directions is investigated. 

2. Nomenclature 

Fx,Fy : Fluid forces in x and y directions 

Fr,Ft : Fluid forces in r and t directions 

Kxx , Ky y , Kxy , Kyx : Stiffness coefficients 

Cxx,Cyy ,Cxy ,Cyx : Damping coefficients 

Mxx,Myy ,Mxy ,Myx : Inertia coefficients 

x,y,z : Fixed coordinates 

r,t : Radial and tangential coordinates 

p(<J>,z) : Pressure distribution 

P(4>) : Average pressure in axial direction 

Pin : Inlet pressure to the seal 

Pex : Outlet pressure from the seal 

AP : Pressure difference between inlet and outlet of the seal 

u : Spinning angular velocity of the rotor 

Q : Whirling angular velocity of the rotor 

e : Whirling eccentricity of the rotor 

$ : Phase difference between principle force and displacement 

R,D : Rotor radius and diameter, respectively 

L : Seal length 

C : Seal clearance 

Vt : Preswirl velocity 

Vts : Preswirl velocity without a rotating motion 

Va : Fluid average velocity in axial direction 

3. Test apparatus and measuring instruments 

3.1 Test apparatus 

Figs. 1 and 2 show the assembly of the test apparatus and the layout of the 
test facility, respectively. In Fig. 1, a working fluid, that is, water, is 

injected through three pairs of swirl passages to accomplish the different inlet 

swirl velocities shown in cross section B. The water passes through the clearances 
of the seal and flows to outlets in both sides of the housing. Cool water is 
continuously applied to the tank in order to maintain constant water temperature. 

The inlet part consists of six tubes connected with swirl passages for 

injecting the water, and swirl speed is adjusted by the six valves in order to 
obtain the arbitrary swirl velocity in the range of 0 - 6.5 m/s. 

The seal assembly consists of a seal stator and a seal rotor. The seal rotor, 
of which the diameter is 70 mm, is connected to the motor by a flexible coupling. 
The motor is controlled in the speed range 0 - 3500 r.p.m. by an electric inverter, 
which also selects the rotational direction. The four types of parallel annular 

seals shown in Table 1 were prepared for the test. The seal stator has holes to 
measure the dynamic pressure in the seal, as shown in cross section C. Both long 
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Table 1 Specifications of seals 


H8M 

Length (mm) 

Diameter (mm) 

L/D 

Clearance (mm) 

qqD 

70 

70 

m 


BiHH 

35 

70 


0.175 

Seal 3 

70 

70 

BO 


Seal 4 

35 

70 

0.5 

0.5 


and short seals have three holes in the axial direction and four holes in the x and 
y directions. The dynamic pressure is measured by the strain gauge type pressure 
gauge shown in cross section C. The fluid dynamic force acting on the stator is 
also directly measured by the load cells shown in cross section C for comparison 
with the data of the pressure transducers. 

The bearing assembly has two ball bearings to make spinning and whirling 
motions. To make a whirling motion, an inside sleeve and an outside sleeve, which 
have a 0.05 mm eccentricity to each other, are attached between the two bearings as 
shown in cross section A of Fig. 1. The two eccentric sleeves can be rotated, 
relatively. So an arbitrary eccentricity can be adjusted in the range of 0 - 0.1 
mm. The sleeves of both sides are driven by a motor through the timing belts. The 
motor can also be controlled by an electric inverter in the rotating speed range of 
0 - 3500 r.p.m., and rotational direction can be selected by an inverter. 

3.2 Measuring instruments 

The measuring procedure illustrated in Fig. 3 consists of four kinds of 
physical variables: that is, rotating and whirling speed, dynamic pressure in the 
seal, seal forces, and displacement of the rotor. Signals from measuring 
instruments are recorded by a data recorder and analyzed by a computer. 

The rotating and whirling speed are measured by eddycurrent type pulse sensors 
and digital counters. 
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A pressure transducer of the semiconductor strain gauge type is used to measure 
the dynamic pressure through the syringe needle. The pressure in the seal is 
measured at the three points in the seal shown in Fig. 4. The seal stator is 
mounted on the housing by 0-rings as shown in Fig. 4. 

Dynamic fluid force can be directly measured by a load cell as shown in Fig. 5. 

The force measured by the load cell is calibrated to revise the influence of the 

0-rings and the inertia of the seal. 

An eddycurrent type displacement sensor is used to measure the displacement, 
and this displacement is used as a reference signal to obtain the phase difference 
between the displacement and the flow induced force. 

The test data are recorded by a data recorder and sent to the computer through 
an A/D converter. In the computer the pressure values are integrated in the 

circumferential direction to obtain the fluid force; then the characteristic 
coefficients are calculated. 

A special pitot tube set is used to measure the preswirl velocity in the seal 
inlet, as shown in Fig. 6. The static and total pressures are measured by the 

pressure transducers instead of the usual U-tube. 

Axial flow velocity in the seal and leakage are determined from the outlet flow 
velocity directly measured by a pitot tube. 

3.3 Calibration 

As a preliminary test, the static and dynamic characteristics of measuring 
instruments (i.e., the pressure transducer, the load cell, and the pitot tube) were 
calibrated . 

Because the pressure measuring set is constructed by a pressure transducer and 
a syringe needle, a phase difference between real pressure and measured pressure was 
investigated. Then, it was found that the pressure measuring set has a phase lag 
relative to the real pressure and that the phase lag increases with the frequency. 
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Fig. 6 Detail of swirl velocity measurement 


Dynamic calibration of the load cell was done in various working conditions. A 
periodic force was excited on the seal by a shaker, and the force was simultaneously 
measured by strain gauges attached to the shaking rod and by the load cell fixed on 
the other side. These data were used to calibrate the measured data in the real 
test. 

A pitot tube set was used to measure the inlet preswirl velocity. Its pitot 
tube coefficient for calibration was determined by means of a standard pitot tube. 


4. Discussions on the accuracy measured by the test apparatus 

4.1 Comparison with the data by two measuring methods 

Fig. 7 illustrates the fluid forces measured by the load cell unit and the 
pressure transducer unit. These results show a reasonable agreement because the 
relative error between the two results is within 5.5 percent. 


4.2 Comparison with other researcher’s results 

A comparison of the authors' experimental results with Childs’ (Ref. 2) results 
is shown in Fig. 8. These results cannot be compared quantitatively, because the 
working fluids and test conditions are different. However, these results have the 
same tendencies. Fig. 9 gives the test results of Kanemori et al. (Ref. 10). These 
results also demonstrate good agreement with the authors' results of Fig. 10. 


5. Calculation of fluid force and characteristic coefficients 

It is assumed that the motion equation of the rotor system is represented as 
follows 


[M]{X}+[C]{X}+[K]{X}={f(t)}-{F} 
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Fig. 8 Comparison of experimental results and Childs results 





Wlilr I/Spin 

Fig. 9 Kanemori's results 



-1.5 -1 -.5 0 .5 1 1.5 

Fr: Ft: - - - Whinl/SpinC 2000r. p. m ) 

Fig. 10 Results for forward and backward spinning speeds 

where [M],[C],[K] are the mass, the damping, and the stiffness matrix and { f ( t ) } is 
the force vector. {F}={Fx Fy} fc is the fluid reaction force acting on a rotor and 
is represented as a linear function of rotor displacement, velocity, and 
acceleration, as follows. 

fFxlfMxx Mxyjfx) fCxx Cxyl(x) fKxx Kxylfx] . . 

(Fyj [MyxMyyj(yJ |cyxCyyJ(y] (KyxKyyj[yJ 

For a small whirling motion about a center position, the relation of the 
coefficients may be expressed by 

Mxx = Myy , Myx = Mxy = 0, Cxx = Cyy, 

Cyx = -Cxy, Kxx = Kyy , Kyx = -Kxy *••(2) 
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where the cross-coupled inertia coefficient is neglected because it is negligibly 
small . 

A rotating coordinate system rotating with the rotor is adopted to analyze the 
fluid forces easily. The relation between a fixed coordinate system and a rotating 
coordinate system is illustrated in Fig. 11, where Fx and Fy are represented by 
the radial force Fr and the tangential force Ft. 



In the test, whirling motion is given by 

x = e«cosflt f y = e*sinfit 
x =-e*fl*sinnt, y = cosftt •••(3) 

x =-e # f2 2 *cosQt, y =-e # ft 2# sinftt 

Substituting Eq. (3) and Eq. (2) into Eq. (1), the following expressions are 
obtained 

Fr = e (-Kxx-ftCxy+ft 2 Mxx) 

Ft = e ( Kxy-fiCxx) 

or 

-Fr/e = Kxx+QCxy-Q 2 Mxx 
-Ft/e =-Kxy+f2Cxx 

where -Fr/e and -Ft/e are called the restitution force coefficient and the 
tangential force coefficient, respectively. 

The displacement of the rotor, the pressure, and the force can be directly 
measured by the displacement sensor, the pressure transducer, and the load cell, 
respectively. So the dynamic coefficients are determined from these data. 

The fluid forces are obtained by integrating the pressure distribution along 
the axial and circumferential directions of the seal: that is, 


•••( 4 ) 

... (5) 
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Fr = /j J**p(<t>>z) •cos<J)*R»d(J)»dz 

= R»L*/”p( 4>) ‘cos^d^ 

= (L*D*7r/N)ZP(<t>j ) »cos<J>j 
) m * 

•••( 6 ) 

Ft = /j'pC'ti 2 ) •sin < J>*R»d<|)»dz 
= R • L • /g*P ( <t> ) • sin<|)«d<|> 

= (L»D»Ti/N).£p(<|>j ) •sin<t>j 


where N is the time series total number of pressures within one period. P($) is 
average pressure in the axial direction. 


P(<t>) = (1/L)/ 0 t p(<j>,z)*dz 
= (1/L||pi»Li 


(7) 


where Pi is the pressure measured at different points, and Li is the length 
between these points. 

To determine the coefficients Mxx, Cxx, Cxy, Kxx, and Kxy of Eq. (5), Fr/e and 
Ft/e at 12 different forward and backward whirling speeds are measured for a given 
spinning speed, and Ft/e is approximated by a linear function, and Fr/e by a 
quadratic function of Q. From Eq. (5), the inertia coefficient, damping 
coefficients, and stiffness coefficients are represented by the curvature, the 
slopes, and the crosses, respectively. The experimental conditions are shown in 
Table 2. 


Table 2 Experimental Conditions 


Pressure differenced P(kPa) 

196,294,490,588,882 

Whirling amplitude e ( u a) 

20,30,40,50,60 

Spinning speed o> (r.p.m) 

500 - 3500 

Whirling speed £2 (r.p.m) 

± 600 ~ ± 2400 

Preswirl velocity Vt (m/s) 

0 - ±13 


6. Experimental results and discussion 

Since the fluid forces measured by the pressure transducer and load cell are 
about the same, the fluid force measured by the pressure transducer is used as the 
experimental results in this discussion. 
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6.1 Pressure distribution and leakage (static results) 

The static pressure distributions in the axial direction of seal 1 and 2 are 
shown in Fig. 12, and the dynamic pressure distributions of seal 1 are shown in Fig. 
13 for various inlet pressures. It is known from Fig. 13 that the dynamic pressure 
distribution is more sensitive to the inlet pressure. 

The leakages of seal 1 to 4 are shown in Fig. 14. They increase as the 
pressure difference and the seal clearance increase but reduce slightly with 
increasing spinning speeds. Also they are insensitive to the preswirls, whirling 
amplitude, and whirling speeds of the rotor. 
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Fig* 12 Static pressure falls 
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Fig. 15 Force coefficients for spinning speed 

6.2 Fluid force 

6.2.1 Effect of spinning speed 

Fig. 15 illustrates the experimental results of fluid force coefficients Fr/e 
and Ft/e for rotor spinning speeds. These results show that Fr/e does not 




obviously vary with the spinning speed, but Ft/e increases with the spinning speed. 
That is, increasing the spinning speed has a destabilizing effect on the rotor 
system. The changes of Fr/e and Ft/e with spinning speed are independent of the 
pressure and the whirling amplitude, and these tendencies appear to be the same for 
all four types of seal. 

The radial and the tangential forces (i.e., Fr/e and Ft/e) for the ratios Q/w 
are shown in Fig. 16 where the oblique fillet is the destabilizing region of Ft/e. 
It is known from the results that Fr/e and Ft/e vary with the spinning speed, even 
if the ratios of fi/co are the same. This may be due to fluid viscosity. 

6.2.2 Effect of preswirl velocity 

Fig. 17 shows the effect of the whirl ratio Q/co on Fr/e and Ft/e of seal 1 for 
various preswirl velocities. The results illustrate that the preswirl in the same 
direction as the spinning affects the unstable area of fluid force greatly. The 
preswirl in the direction of rotation extends the unstable area, but the preswirl in 
the direction opposite to that of rotation reduces the area, even reducing it to 
negative values for enough opposing swirl. However, for both of the swirls the 
intersection point of Ft/e and the zero axis approach about Q/co = 0.5 with 
increasing spinning speed. 

This phenomenon is illustrated as the relation between the whirling speed and 
the average circumferential velocity rotating about coR/2. That is, when the 
whirling angular velocity of the rotor is lower than the average angular velocity of 
the fluid, the whirling motion of the rotor will be accelerated by the rotating 
fluid; conversely, when the whirling angular velocity is faster than the average 
angular velocity, the whirling motion will be restrained, so that the sign reverses 
about Q/co = 0.5. The circumferential velocity distributions in radial are shown in 
Fig. 18 where Figs, (b) and (c) are for the cases of enough preswirl and inverse 
preswirl velocities. The distribution demonstrates that the preswirl strengthens 
the fluid average velocity Vm, and the opposing preswirl weakens the average 
velocity Vm. From these effects, the unstable area varies. 

In summary, the preswirl in the direction of the rotor rotation has a 
destabilizing effect on the rotor system, and the opposing preswirl has a 
stabilizing effect on the rotor system. 

6.2.3 Effect of pressure difference 

Fr/e and Ft/e are measured for various pressure differences between the inlet 
and outlet of the seal. The results show that Fr/e and Ft/e only change their 
values; they do not change their qualitative characteristics. One of these results 
is shown in Fig. 19. 

6.2.4 Effect of seal length 

The effect of seal length is shown in Fig. 20. The result shows that the two 
seals have the same tendencies for Ft/e, but not for Fr/e. Compared with that of 
the long seal (L/D = 1), the Fr of the short seal (L/D = 0.5) acts as a stabilizing 
force. This is because if the axial flow velocity increases, the tangential force 
which is generated by the circumferential flow decreases with the decrease of the 
time it takes for the flow to pass through the seal. 

6.2.5 Effect of whirling amplitude 

Fig. 21 shows the effect of whirling amplitude (eccentric ratio) on Fr/e and 
Ft/e of seal 1. The results show Fr/e and Ft/e increase with increasing whirling 
amplitude. 


242 



Ft/e (kN/mm) . . Fr/e (kN/mm) 


5 

2.5 
0 

2.5 
5 

-2 -1.5-1 -.5 0 .5 1 1.5 2 

Whirl/Spin 



** Spln-1500r. p. n. 
++ Splrv*2500r. p. m. 
oo Spln-3500r. p. m. 



AP=882 kPa 
Vts=0 m/s 
e=0.05 mra 

i 

1 



(a) Radial coefficient Fr/e 



-2 -1.5-1 -.5 0 .5 1 1.5 2 

Whirl/Spin 


(b) Tangential coefficient Ft/e 
Fig. 16 Effect of spinning speed on force coefficients 
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Fig. 17 Effect of swirling velocity on force coefficients 


244 






Stator 

l 


Vm 




r~ 



Rotor 

| coR 

1 


(b) Positive swirl 


Stator 

— 

Vm 




Rotor 

<uR 



(c) Hegative swirl 


Fig. 18 Velocity distributions in clearance of seal 
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Fig. 19 Effect of pressure difference on force coefficients 


6.2.6 Effect of seal clearance 

The results of Fr/e and Ft/e for different radial seal clearances are shown in 
Fig. 22. They are measured in the condition of positive preswirl velocities. 
Attention should be paid to Fr/e; that is, regardless of the seal length (L/D = 1.0 
and L/D = 0.5), seals with large clearance have negative values of Fr/e, in contrast 
to seals with small clearance. It is considered that these negative Fr/e are caused 
by the increase of the axial flow in the large clearance seal. In addition, the 
tangential force coefficients Ft/e of the large clearance seals are much more 
influenced by preswirl velocity than those of the small clearance seal. 
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6.3 Stiffness, damping, and inertia coefficients 

The stiffness coefficients, damping coefficients, and inertia coefficients 
obtained from the fluid forces are shown in Figs. 23 and 24. These results show 
that the inertia coefficients are independent of the pressure difference, the 
eccentric ratio, and the swirl velocity, and that they always have positive values. 
Therefore, the inertia force decreases an eigenf requency of the rotor system. 

The direct damping coefficients are always positive and much larger than Cxy, 
so generally, the damping forces act as stabilizing forces on the rotor system. 
These damping forces are not affected by swirl velocity and seal length. 

The stiffness coefficients Kxx and Kxy are largely affected by the swirl 
velocity, especially Kxy. These can be recognized from Fig. 23. Kxy increases with 
the preswirl velocity and decreases with the opposing swirl velocity. That is, the 
stiffness force related to Kxy may act as a stabilizing force or as a destabilizing 
force according to opposing preswirl or preswirl velocities. 

As shown in Fig. 24, the coefficient Kxx has a negative value for the long 
seal 1, but Kxx has a positive value for the short seal 2. Namely, the short seal 
is more stable than the long seal from the view point of the stiffness coefficients 
of diagonal terms. 


7. Conclusions 

From the experimental study reported here, we derive the following conclusions: 

(1) Leakage grows with increasing pressure difference, but reduces slightly with 
increasing spinning speed. The leakage is insensitive to the preswirl velocity, the 
whirling amplitude, and the whirling speed. 

(2) Flow induced force increases with increasing whirling amplitude, but decreases 
with increasing seal clearance. 

(3) Preswirl velocity strongly affects the stability of the rotor system. Preswirl 
velocity in the rotating direction has a tendency to destabilize, but opposing 
preswirl has a tendency to stabilize the rotor system. 

(4) Tangential fluid force increases with rotor spinning speed; as a result, it may 
become a destabilizing force. 

(5) Axial flow through the seal has a tendency to stablize the rotor system. It can 
be recognized from the radial force coefficients Fr/e of the short seal and the 
large clearance seal. 

(6) The inertia force related to the inertia coefficient Mxx always acts to 
decrease the eigenf requency of the rotor system. 

(7) The damping force related to the damping coefficients Cxx and Cxy usually acts 
as a stabilizing force. As for the values, Cxx is larger. 

(8) The preswirl velocity has no effect on Mxx, Cxx, or Cxy; however, it affects the 
stiffness coefficients Kxx and Kxy greatly, especially Kxy. 

The financial support of Grant-in-Aid for Scientific Research B, Ministry of 
Education, and also Aero-Engine and Space Operation, Ishikawajima-Harima Heavy 
Industries Co. Ltd. is gratefully acknowledged. The authors wish further to thank 
Dr. S. Saito for his kind suggestions for this research. 
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